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Abstract 

We present the result for the two-loop and the one-loop squared virtual QCD corrections to the 
W boson pair production in the quark-anti-quark-annihilation channel in the limit where all kine- 
matical invariants are large compared to the mass of the W boson. The infrared pole structure is 
in agreement with the prediction of Catani's general formalism for the singularities of two loop 
amplitudes. 



1 Introduction 



The Large Hadron Collider (LHC) will be the centre of interest for particle physics phenomenology 
in the next years. Open issues that require definite answers are the verification of the consistency 
and validity of the Standard Model (SM) in the energy range of the LHC as well as insights into 
New Physics. Several proposed models and concepts that have the SM as their low energy limit 
theory are either to pass the LHC test or to be proven wrong. Supersymmetry and Extra-dimensions 
are two of the most illustrious examples. 

Probably, the most important goal for the LHC is the discovery of the elusive Higgs boson. 
The latter is part of the mechanism of dynamical breaking of the Electroweak (EW) symmetry and 
is responsible for the fermions and gauge bosons mass. Discovering the only constituent of the 
Standard Model (SM) which has not been experimentally observed yet, along with a systematic 
measurement of its properties, will be essential for our understanding of mass and the precise gauge 
structure of the SM. Another important endeavour at the LHC, in connection to the investigation of 
the non-Abelian gauge structure of the SM, is the precise measurement of the hadronic production 
of gauge boson pairs, WW, WZ, ZZ, Wy, Zy. Deviations from the SM predictions would indicate 
the presence of either anomalous couplings or new heavy particles which would decay into vector 
boson pairs [1,2]. 

Seen under the prism of the previous argumentation, W pair production via quark-anti-quark- 
annihilation, 

qq^>W + W~, (1) 

is a very important process at the LHC. Firstly, it can serve as a signal process in the search for 
New Physics since it can be used to measure the vector boson trilinear couplings as predicted 
by the Standard Model (SM) (actually, this is the favored channel as it involves both trilinear 
vertices, WWZ and WWy). Secondly, qq — > is the dominant irreducible background to the 

promising Higgs discovery channel 

pp^H^ W*W* -> IvFV , (2) 

in the mass range MHi ggs between 140 and 180 GeV [3]. 

Due to its importance, the study of W pair production in hadronic collisions has attracted a lot of 
attention in the literature. The Born cross section was calculated almost 30 years ago [4], whereas 
the next-to-leading order (NLO) QCD corrections to the tree-level were computed in Refs. [5-9] 
and were proven to be large. They enhance the tree-level by almost 70% which falls to a (still) 
large 30% after imposing a jet veto. Therefore, if a theoretical estimate for the W pair production 
is to be compared against experimental measurements at the LHC, one is bound to go one order 
higher in the perturbative expansion, namely to the next-to-next-to-leading order (NNLO). This 
would allow, in principle, an accuracy of better than 10%. Notice that first steps in this direction 
have been done by considering soft-gluon resummation effects in W pair production [10]. 

High accuracy for the W pair production is also needed when the process is studied as back- 
ground to Higgs production. The NLO QCD corrections to the signal process for the Higgs dis- 
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covery via gluon fusion, gg — > H, contribute a 70% [11, 12], whereas the NNLO contributions 
suggest an additional 20% for the LHC [13-15]. With a jet veto, at NNLO the total corrections 
are of the order of 85% [16-18]. Lastly, the QCD corrections to the cross section for the pro- 
cess H — > WW — > l\Tv' are known at NNLO [19, 20] whereas the EW ones are known beyond 
NLO [21]. The ratio of the Higgs signal over background is expected between 1:1 and 2:1 once 
certain cuts are applied that reject events with high pr jets. For a consistent QCD analysis, there- 
fore, we need to compare both signal and background cross sections calculated at the same order, 
that is, at NNLO. Another process that needs to be included in the background is the W pair pro- 
duction in the loop induced gluon fusion channel, 

gg^W+W-. (3) 

This contributes at O(oc^) relative to the quark-anti-quark-annihilation channel but is nevertheless 
enhanced due to the large gluon flux at the LHC. The corrections from gluon fusion increase 
the W pair background estimate by almost 30% after certain experimental Higgs search cuts are 
imposed [22,23]. 

In this paper, we address the task of computing the NNLO virtual part, more precisely the inter- 
ference of the two-loop with the Born amplitude, as well as the the one-loop squared contribution. 
We work in the limit of fixed scattering angle and high energy, where all kinematical invariants 
are large compared to the mass m of the W. Our result contains all logarithms logm as well as 
all constant contributions while we neglect power corrections in m. These will be presented in a 
following publication. 

Our methodology for obtaining the massive amplitude (massless fermion-boson scattering was 
studied in Ref. [24]) is very similar to the one followed in Refs. [25-27] which is, at its turn, an 
evolution of the methods employed in Refs. [28,29]. The amplitude is reduced to an expression 
that only contains a small number of integrals (master integrals) with the help of the Laporta al- 
gorithm [30]. In the calculation for the two-loop amplitude there are 71 master integrals. For the 
one-loop squared case, we use the helicity matrix formalism to reduce the problem to a small set of 
integrals. Next comes the construction, in a fully automatised way, of the Mellin-Barnes (MB) rep- 
resentations [31,32] of all the master integrals by using the MBrepresentation package [33]. The 
representations are then analytically continued in the number of space-time dimensions by means 
of the MB package [34], thus revealing the full singularity structure. An asymptotic expansion in 
the mass parameter is performed by closing contours and the integrals are finally resummed, either 
with the help of XSummer [35] or the PSLQ algorithm [36]. 

Our paper is organised as follows. In Section [2] we introduce our notation, present briefly our 
methods and define the perturbative expansion of the matrix elements summed over colours and 
spins. In Section [3] we study the singular behavior of the NNLO contributions, and verify that 
it agrees with the general formalism developed by Catani [37] for the infrared structure of QCD 
amplitudes. In Section |4] we present the finite remainders for the interference of the tree and the 
two-loop amplitude and the one-loop squared after subtraction of the singular poles of Section [3] 
from the explicit result. We organise the finite part according to the colour content of the two-loop 
amplitude for the two-loop case. The finite remainders are expressed in terms of logarithms and 
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polylogarithms which are real in the physical domain. We conclude in Section [5l Finally, for 
completeness, the one-loop result up to order £ 2 is included in the Appendix. 



2 Notation 



The charged vector-boson production in the leading partonic scattering process corresponds to 

qj(pi) + qj(p2) -> W-(p 3 ,m)+W + (p4,m), (4) 

where p, denote the quark and W momenta, m is the mass of the W boson and j is a flavour index. 
We are considering down type quark scattering in our paper. Obtaining the corresponding result 
for up-type quark scattering is actually trivial as we will show in the following. Energy-momentum 
conservation implies 

P" + P2=P3+P4- (5) 

We consider the scattering amplitude M for the process © at fixed values of the external parton 
momenta pu thus p\ = p\ = and p 2 = p\ = m 2 . The amplitude M may be written as a series 
expansion in the strong coupling a s , 
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M) = h (0))+ (S)^ (1))+ (S) ^ (2) >+°^) 



(6) 



and we define the expansion parameter in powers of 0C S (// 2 ) / (2jt) with /j being the renormalisation 
scale. We work in conventional dimensional regularisation, d = 4 — 2e, in the MS-scheme for the 
coupling constant renormalisation. 

We explicitly relate the bare (unrenormalised) coupling to the renormalised coupling a s by 



(7) 



where we set the factor S £ = (47l) e exp(— eje) = 1 for simplicity and [3 is the QCD [3-function 
known at present up to the four-loop level [38, 39] 

Po = yQ-^7>n / . (8) 

The color factors in a non-Abelian SU(iV) -gauge theory are C A =N,C F = (N 2 — 1)/2N and T F = 
1/2. Throughout this paper, Af denotes the number of colors and n f the total number of flavors 
of massless quarks. Remark, however, that the latter must come in pairs, because of the flavor 
changing coupling to the charged gauge boson. This is only problematic in the case of top quarks 
running in a closed loop. 

In the following, our discussion will be restricted to the two-loop amplitude summed over 
spins and colours and contracted with the Born one. Nevertheless, it should be stressed that our 
methods and the results of the present work can be easily extended to the partial amplitudes for the 
individual helicity combinations of the massive two-loop amplitude \ M^) itself. 
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For convenience, we define the function A (e, m, s, t , /j) for the squared amplitudes summed over 
spins and colors as 



Y,\M(qj + qj^W++W-)\ 2 = A(E,m,s,t,n). 
SI is a function of the Mandelstam variables s, t and u given by 

ry ry ry ry ry 

s={p\+P2) , t = (pi -pi) -m , u = (pi-p 4 ) 
and has a perturbative expansion similar to Eq. ©, 



m 



Ji(e,m,s,t,/j) 



In terms of the amplitudes the expansion coefficients in Eq. (fTTI) may be expressed as 

,(0) = (^(0)^(0)^ 



A 1 



(2) 



(9) 
(10) 

(11) 

(12) 
(13) 

(14) 



where M (°) and fW (*) are the massive tree level and one loop amplitudes correspondingly. J? is 
given by 



(0) 
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+C 3 
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( 1 — x) m s m L s 



,„ \ /1 (3-4e)-2x(l-x) 4x(l-x) 
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m. 



xx -. ^ /„ /- xx 6 — 8£ — 8x(l — x) 2x(l—x) 
-24(1 -x(l -jc)) + 16e (2— jc(1 —x)) H H v „ ; 



m. 



(15) 



where we have defined x = m s = IJ j- and only the leading physical powers (i.e. down to 
the constant) in the m s -expansion are retained. Notice that, once the actual values of the c, are 
substituted, the terms singular in m s cancel as required by unitarity. This will be the case for the 
final two-loop and one-loop squared expressions as well. The coefficients c\, C2 and C3 are in their 
essence combinations of EW coupling constants defined as 



SWL 



In order to compute the (M^\M^>) it proves convenient to express the amplitude \M^) 
in terms of helicity amplitudes, M g {X\,X 2 ,s,t), where X\ and X 2 stand for the helicities of the 
W + and W respectively. In other words, it is convinient to decompose the amplitude into a 
sum of products consisting generally of three parts, a Feynman integral, a rational function of the 
kinematical variables and a standard matrix element, M a complete list of which is listed below in 
Eq. (fT8l) . For that, Dirac algebra is used, as well as the equations of motion and an anticommuting 
Y5. The quark and anti-quark have opposite helicities in the centre-of-mass system so one helicity 
label above, g = ±1, suffices. 

Therefore, the one-loop amplitude is formally rearranged as 

= ^QisS^l/isJiU-^Mjdpk}^), (17) 

where the Q are coefficients, the if are one-loop dimensionally regularized scalar integrals, Mj are 
helicity matrix elements, g = ± and k = 1, ...,4. The ten helicity matrix elements Mj(pk,g) = M? 
have been taken as defined in Ref. [40] (see also [41]): 
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= v{ P2 )i 


\{p 3 -]/> 2 )t 2 'p g <p 




Ml 


= v(p 2 )f 


yp g w(>i)ei -£ 2 , 




Ml 


= v(pi)i 


l ( P g u{pi)Z 2 -p3, 




Ml 


= -V(p 2 , 


)^gu(pi)e 1 -p 4 , 




Ml 


= v(Pi)t 


l Vgu(pi)e 2 -pi , 




Ml 


= ~V(P2, 


)^gU(pi)E 1 -p 2 , 




Ml 


= v(p 2 )$ 


< 3 P g u(pi)£i-p 2 £ 2 - 


Pi, 


Ml 


= v(p 2 )$ 


< 3 P g u(pi)Ei -p 2 e 2 - 


P3, 


Ml 


= v(p 2 )f 


> 3 <P g u(pi)Ei -P4Z 2 - 


Pi, 


Ml 


= v(p 2 )f 


< 3 V g u(pi)£i -P4E 2 - 


P3, 



where <P g = T± = All colour indices as well as the arguments of the polarization vectors, 

£i(/?3,Xi) and £2(p4,A,2), have been suppressed. 

After expressing the one-loop amplitude as in Eq. (fTTb and calculating the Feynman integrals, 
it is trivial to obtain (M ^>\M W), one needs only to compute the traces (one fermionic chain 
in all cases) coming from multipling a matrix element with the complex congugate of another 
one. We have decomposed the tree level amplitude as well in terms of helicity amplitudes and 
computed (M ^ \M ^) as a trivial cross check. Even though the representations in Eq. (TT8l) have 
been used internally, here we present our result only for the amplitude squared and summed over 
helicities. Notice that this is done in conventional dimensional regularization, which implies 2 — 2£ 
polarizations of the vector bosons. 

The expressions for & V-> have been presented e.g. in Refs. [5, 6] whereas the leading color 
coefficient of (M ^ \M ^) was discussed in Ref. [42]. Here we provide for the first time the result 
for the real part of R ^ . 
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3 Infrared Pole Structure 



In the simpler case of one-loop amplitudes, their poles in e can be expressed as a universal combi- 
nation of the tree amplitude and a colour-charge operator ^^(e). The generic form of the ^(e) 
operator was found by Catani and Seymour [45] (see also [43, 44]) and it was derived for the 
general one-loop QCD amplitude by integrating the real radiation graphs of the same order in 
perturbation series in the one-particle unresolved limit. 

The pole structure of our one-loop expression is given, according to the prediction by Catani, 
by acting with the operator I^ 1 ) (e) onto the tree-level result: 

\af®) = I (1) (e)l^ (0) ) + I^ fi ( l), (19) 

where I*-^(e) is defined as 



r(l-e) Ve 2 2e 

In a similar way, the divergences of the two-loop amplitude can be written as a sum of two 
terms: the action of the I^ 1 ) (e) operator on the one-loop amplitude and the action of a new operator 
I^ 2 ) (e) on the tree amplitude. The I^ 2 ) (e) operator includes a renormalisation scheme dependent 
term H^> multiplied by a 1/e pole. In the following, we give explicit expressions for T^(e) and 
l( 2 )(e) which are valid in the MS scheme. 

At next-to-next-to-leading-order (NNLO), contributions from the self-interference of the one- 
loop amplitude and the interference of the tree and the two-loop amplitude must be taken into 
account, so that 

A^(j,M,m,Ai) = A^< lxl >(j^^ (21) 

^ NNL °( lxl )(5,f, M ,m, i u) = (^( 1 )|5W-W), (22) 
^NNLO(0x2) (j)?)M)m ^ )= (^(0)1^(2)^^(2)^(0)^ (23) 



with 
and 



We further decompose the one-loop self-interference and the two-loop contributions as a sum 
of singular and finite terms, 

^NNLO(lxl) (S;M)m;iU) = C j}£\s,t,u,m,ti) + !F&?\s,t,u,m,ti) (24) 

and 

*NNIX)(0x2) (j|M)OT ^ ) = c M (25 ) 



cjttani an d Catani contam infrared singularities that will be analytically canceled by the in- 
frared singularities occurring in radiative processes of the same order (ultraviolet divergences hav- 
ing already been removed by renormalisation). ^J^* and are the remainders which are 
finite as £ — > 0. 
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The infrared poles of the interference of the tree and the two-loop amplitudes follow a generic 
formula developed by Catani in Ref. [37]. Due to the simple colour structure of the process © the 
action of lW(e) and l( 2 )(£) is factorised such that we formally have 

i^M^Hl^ (26) 
dZf fa, u, m^) = 2Re {iW (e) (M® |M«) + 1^ (e) (M® \M® > } (27) 



and 



with 



**>(«) = 4 I «..(e)(l<..(e) + f) + £^p(|«)l<..( 2 e) 

+// (2) (£), (28) 



where 



67 7C 2 \ „ 10. 



K= iw-T CA -T T '"f <29) 



The renormalisation scheme dependent constant for a QCD amplitude with a qq pair is given 
by 

^+£-3*)<M-l + *W- <3o) 

We were able to verify that our results have the same infrared structure as the one predicted by 
Catani's formalism. 



4 Results 

4.1 Two-loop Contribution 

In this section, we give explicit expressions for the finite remainder of the two-loop contribution 
finite defined as 



^fV^^H*^ ^^ (3D 

or in the rescaled form 
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s ■) 


- 7 (° X2 ) ( m y ±) 
H 




s \ 

M 


(0x2) , S 

~ h,down\ m ^y-> 2 
A* 




s \ 


— J 3,down\ m ^y^ 2> 
H 




s 


I J A,down\ m ^y^ 2' 



The EW structure of the finite remainder for a down-type quark can be factorised as 

= 2N I Ci&rn*^- (33) 

(=1,4 r 

This decomposition allows one to easily obtain the result for the up-type quark scattering. The 
latter is then given by 

f£l = 2iV E c i3 Vf\m„x°), (34) 

1=1,4 t* 

where one needs to use the following formulae 

(35) 



(36) 
(37) 
(38) 

and naturally to make the corresponding changes in the definitions of the couplings c\, C2, C3 and 
C4, namely to use the up-type quark charge and isospin. Here y = —-. In the following and with 
no loss of clarity, since our result assumes down-type quark scattering, we will suppress all indices 
that indicate the type of scattered quark. The functions Ji(m s ,x, A) in Eq. (1331) will be presented 
decomposed according to the colour structure, namely in the form 

4 0x2 \m s ,x,^- 2 ) = (jl l k F C A + j^ F + jPc F T F n f ). (39) 

C4, in addition to c\, C2 and C3, is a new coupling that appears at the two-loop level and is defined 

as: 

= C«g"zA {m 

24(1 -f) 

The appearance of C4 is an effect that comes from a specific part of (M ^>). This part 
consists of two-loop fermionic boxes contracted with the Born diagram that involves an ^-channel 
Z exchange. A typical example can be seen in Fig. (1). The main feature of these diagrams is that 
their EW couplings fall into two disjoint fermionic chains and once the traces are computed the 
axial part drops out. By adding and subtracting to the surviving vector part the corresponding axial 
part, one can combine vector and axial contributions into a piece proportional to 02- The remaining 
piece is proportional to what we have defined as C4. 
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Fig 1 : Born diagram with a Z exchanged in the s-channel contracted with a fermionic two-loop box. 

We have verified that applying the naive recipe of sending all traces that contain a single 75 
independently to zero is a valid approach for this class of diagrams. We did this by calculating ex- 
plicitly the output after substituting 75 by its alternative form 75 = ^T£ it/ vapY^Yv7a7p an d confirming 
that no additional terms survive. This was in fact a non-trivial cancellation as it occurs only for the 
sum of all the diagrams of this particular class. Note, however, that there were finite contributions 
from traces containing 75 in the case of pure W boson pair exchange (no photons or Z's involved). 

We are finally ready to present our result. The functions 3u are given by 
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where L m , L., Ly and L v are defined as 



L m = log (hi,) , L 5 = log (j£ j , L x = log (x) , L y = log (1 - x) . 



(45) 



4.2 One-loop Squared Contribution 

In this section, we give explicit expressions for the finite remainder of the one-loop squared con- 
tribution ^ defined as 



(46) 
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The EW structure of the finite remainder for the one-loop squared corrections, similarly to the 
case of the two-loop corrections, can be factorised as 
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Our result then reads: 
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5 Conclusions 

In this work we have calculated the NNLO QCD virtual corrections for the process qq — > W + W~ 
in the limit of small vector boson mass. The MS renormalised amplitude is still infrared divergent 
and contains poles up to 0(1/e 4 ). We checked that the infrared structure of our result agrees with 
the prediction of Catani's formalism for the infrared structure of QCD amplitudes. 

The main result of our paper has been given as the finite remainder of the NNLO two-loop 
and one-loop virtual corrections after subtraction of the structure predicted by Catani's formalism. 
This is a first step towards the complete evaluation of the virtual corrections. In a forthcoming 
publication, we will derive a series expansion in the mass and integrate the result numerically. 
This will require the present result as a starting point. 

To complete the NNLO project one still needs to consider 2 — > 3 real-virtual contributions 
and 2 — > 4 real ones. The real- virtual corrections are known from the NLO studies on WW + jet 
production in Refs. [46,47]. The integration over the full phase space would require additional 
subtraction terms, similar to those constructed in Ref. [48]. 
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Appendix: (M (°) | M^^) to order e 2 

Here we present the expression for the one-loop result, (M ^|^W^j te ) up to order £ 2 for down- 
type quarks. This result completes the list of the elements needed in Eq. (fTTI) in order to have the 
perturbative expansion of the amplitude up to order o? s in the high energy limit. 
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-8 (-32x 2 + 32x+ (5x 2 -5x+ l) & - 16) + 16 (-8x 2 + 8x+ (3x 2 - 3x + 3) £3-4) L,] } 
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